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ABSTRACT 

We explore the nature of unresolved X-ray emission in a broad sample of galaxies of all mor- 
phological types based on archival Chandra data. After removing bright compact sources, 
we study LxILk luminosity ratios of unresolved emission, and compare them with the So- 
lar neighborhood values. We conclude that unresolved emission is determined by four main 
components, three of which were known before: (i) The population of faint unresolved sources 
associated with old stellar population. In early-type galaxies, their 2-10 keV band luminosity 
scales with the stellar- mass with LxILk - (3.1 + 0.9) x 10^^ erg s"' Lj^^; (ii) The ISM with 
kT ~ 0.2 - 0.8 keV present in galaxies of all types. Because of the large dispersion in the 
gas content of galaxies, the size of our sample is insufficient to obtain reliable scaling law 
for this component; (iii) The population of unresolved young stars and young stellar objects 
in late-type galaxies. Their 2-10 keV band luminosity scales with the star-formation rate 
with Lx/SFR ^(1.1 + 0.9) x 10^** (erg/s)/(M0/yr); (iv) In four old and massive Virgo ellip- 
ticals (M49, M60, M84, NGC4636) we find anomalously high X-ray emission in the 2 - 10 
keV band. Its presence has not been recognized before and its nature is unclear. Although it 
appears to be stronger in galaxies having stronger ISM component, its existence cannot be 
explained in terms of an extrapolation of the warm ISM spectrum. Association with Virgo 
cluster of galaxies suggests that the excess emission may be due to intracluster gas accreted 
in the gravitational well of a massive galaxy. We investigate this and other possibilities. 

Key words: Galaxies: elliptical and lenticular, cD - Galaxies: irregular - Galaxies: spiral - 
X-rays: galaxies - X-rays: stars - X-rays: ISM 



1 INTRODUCTION 

The X-ray appearance of most nearby (D < 30 Mpc ) galax- 
ies is determined by bright X-ray binaries jpabbianol 12006) . such 
as high-mass X-ray binaries (HMXBs) and low-mass X-ray bi- 
naries (LMXBs). Because of their high X-ray luminosity (Lj ~ 
10^^ - 10^' erg s"') notable fraction of them can be resolved with 
moderately deep Chandra exposur es, and hence hence they can 
be studied in full partic ulars (e.g. iGrimm et alj l2003l : iGilfanovl 
l2004lMineo etalj|201ll) . Underneath bright X-ray binaries dif- 
fuse emission is present in all typ es of galaxies. S imilarly to the 
Galactic Ridge X-ray emission ( Revnivtsev et al. 

I [2006), a part 

of this emission originates from the collective emission of faint 
{Lx ~ 10-^ - 10^^ ergs"') compact X-ray sources, mostly from 
active binaries (ABs), cataclysmic variables (CVs), young stars, 
and young stellar objects (YSOs). Based on the analysis of a few 
nearby early-type galaxies, it was suggested that the overall lumi- 
nosity of ABs and CVs is tightly correlated with the stellar mass of 
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the galaxy (e.g. ^ Pellegrini & Fabbiano|[l994l : iBogdan & Gilfanovl 
120081 : iRevnivtsev et al.ii2008l) . Additionally, in all types of galaxies 
truly diffuse emission from warm {kT < 1 keV) ionized interstel- 
lar medium (ISM) is also present, whose measure, among others, 
depe nds on the mass of the host galaxy, albeit with larg e disper- 
sion ( IO'SullivaneTai]|200ll : iMathews & BrighentillloOa . The to- 
tal observed unresolved emission is the combination of these, and 
possibly other, components. 

Although during the operation of Chandra large number of 
nearby galaxies were observed, a systematic study of the unre- 
solved X-ray emission we re only per formed in a few of the m (e.g. 
iBogdan & Gilfanovl |2003) . Recently, iBoroson et all ( 1201 ih inves- 
tigated a sample of early-type galaxies, mainly focusing on their 
LMXB population and ISM content. The main motivation of the 
present paper is to study the unresolved emission in a large sample 
of early-type and late-type galaxies in a uniform manner, thereby 
gaining a better insight into the importance of various unresolved 
X-ray emitting components. 

Our primary goal is to measure X-ray-to-K-band luminosity 
ratios {LxILk) of the unresolved emission in the 0.5 - 2 keV and 
in the 2-10 keV energy ranges and compare their values across 
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Table 1. The list of early-type and late-type galaxies studied in this paper. 



Name 


Distance 


Lk 


Nh 


Morphological 


SFR 


Age 




Tflit 


^lim 


R 




(Mpc) 


(Lk.o) 


(cm-2) 


type 


(Moyr-') 


(Gyrs) 


(ks) 


(ks) 


(erg s"') 


(") 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


M51 


8.0" 


6.4 X lO'" 


1.6 X 10^" 


SAbc 


3.9 




90.9 


80.8 


9 X 10^* 


175 


M74 


7.3" 


1.6x lO'" 


4.8 X 10^" 


SA(s)c 


1.1 




104.4 


84.7 


7 X 10^* 


140 


M81 


3.6* 


4.8 X 10'" 


4.1 X 10^" 


SA(s)ab 


0.4 




239.1 


202.1 


8 X 10'^ 


250 


M83 


4.5" 


3.9 X 10'" 


3.9 X 10^" 


SAB(s)c 


2.8 




61.6 


52.1 


4 X 10^* 


200 


M94 


4.7" 


3.3 X 10'" 


1.4 X 10^" 


(R)SA(r)ab 


1.2 




76.9 


69.5 


4 X 10^* 


164 


M95 


10.1" 


3.3 X 10'" 


2.9 X 10^" 


SB(r)b 


1.1 




120.1 


99.6 


1 X 10^' 


88 


MlOl 


7.4" 


3.4 X 10'" 


1.2 X 10^" 


SAB(rs)cd 


1.3 




1070.6 


833.3 


6 X 10'^ 


167 


NGC2403 


3.2* 


5.0 X lO' 


4.2 X 10^" 


SAB(s)cd 


0.4 




224.0 


184.2 


7 X 10^^ 


140 


NGC3077 


3.8" 


2.6 X lO' 


4.0 X 10^" 


10 pec 


0.3 




54.1 


42.1 


3 X 10^^ 


77.5 


NGC4214 


2.9" 


6.5 X 10* 


1.5 X 10^" 


lAB(s)m 


0.2 




83.4 


53.6 


1 X 10'^ 


70 


NGC4449 


4.2" 


3.4 X 10' 


1.4 X 10^" 


IBm 


0.4 




102.1 


97.0 


2 X kP^ 


95 


M31 bulge 


0.78" 


3.7 X 10'" 


6.7 X 10^" 


SA(s)b 




5.1^ 


180.1 


143.7 


2 X kP' 


720 


M32 


0.805* 


8.5 X 10** 


6.3 X 10^" 


cE2 




3.8-'' 


178.8 


172.2 


1 X kP* 


90 


M49 


16.3'' 


3.0 X lO" 


1.7 X 10^" 


E2 




8.5-^ 


60.9 


55.7 


5 X 10^' 


140 


M60 


16.8'" 


2.2 X lO" 


2.2 X 10^" 


E2 




11. 0'^ 


109.3 


90.5 


2 X 10^' 


125 


M84 


18.4'' 


1.6 X lO" 


2.6 X 10^" 


El 




11.8^ 


117.0 


113.8 


3 X 10^' 


100 


M89 


15. S'' 


7.0 X lO'" 


2.6 X 10^° 


EO- 




9.6^^ 


55.1 


51.8 


4 X 10^' 


80 


M105 


9.8^* 


4.1 X lO'" 


2.8 X 10^" 


El 




9.3' 


341.4 


314.0 


2 X 10'* 


90 


NGC1291 


8.9'' 


6.3 X 10'" 


2.1 X 102" 


(R)SB(s)0/a 






76.7 


51.2 


1 X 10" 


130 


NGC3377 


11.2'' 


2.0 X 10'" 


2.9 X 102" 


E5-6 




4.1/ 


40.2 


34.0 


3 X 10" 


79 


NGC3585 


20.0'^ 


1.5 X 10" 


5.6 X 10^" 


E7/S0 




3.1/ 


95.9 


89.1 


3 X 10" 


180 


NGC4278 


i6.r 


5.5 X 10'" 


1.8 X 10^" 


El-2 




10.7/ 


467.7 


443.0 


4 X lO^* 


110 


NGC4365 


20.4'' 


1.1 X 10" 


1.6 X 10^" 


E3 




3.6S 


198.3 


181.4 


2 X 10" 


80 


NGC4526 


16.9'' 


9.6 X 10'" 


1.7 X 10^" 


SAB(s)0" 




1.7* 


44.1 


34.5 


5 X 10" 


84 


NGC4636 


14.7'' 


8.1 X 10'" 


1.8 X 10^" 


E/S0_1 




10.3' 


212.5 


202.1 


4 X 10" 


102 


NGC4697 


11.8'' 


5.1 X 10'" 


2.1 X 10^" 


E6 




8.2/ 


195.6 


162.0 


8 X lO'* 


95 



Note. Columns are as foll ows. (1 ) Refer enc es are: "iKarachentsev et al.' ('2004") - * Freedman et al feOOlS - '^[Leonard et alj j2004h - °[ Stanek & Gamavict] 
(1998); Macri ( 1983) - * lMated j 19981) - ^ iTonrv e t al. 1 2001 ) - Jensen et al. 1 2003) - de Vaucouleurs j 19751) . (2) Total near-infrared luminosity of the 
elliptic region described in column (8). (3) Galactic absorption (Dickev & Lockmanlll99 Q). (4) Taken from NED (http://nedwww.ipac.caltech.edu/ I. (5) 
Star-formation rate for late-ty pe galaxies within the same elli pse. (6) Age of the ste llar population for early-ty pe galaxies. For NGC1291 no age reference has 
been found. References are: / |Terlevich & Forbed <2002h - s benicolo et^ j2005h - * ballagher etaP <2008l) - ' ISanchez-Blazguez et al . ( 2006). (7) and (8) 
Exposure times before and after flare filtering. (9) Source detection sensitivity in the 0.5 - 8 keV energy range. (10) Major axis of the studied elliptic regions. 
The orientation and shape of the regions were taken from K-band measurements (http://irsa.ipac.caltech.edu/applications/2MASS/). 



morphological types and stellar masses. In particular, we aim to 
investigate whether Lx/^a: ratios of unresolved compact sources 
in the old stel lar population are in deed universal as was suggested 
earlier by Rev nivtsev et al ] (l2008h . Extending this investigation to 
late-type galaxies we intend to study LxlLic ratios of a stellar pop- 
ulation that also includes young stars and YSOs. Thus, the X-ray 
emission from these sources and their contribution as a function 
of the star-formation rate (SFR) can be studied. Furthermore our 
sample includes low-mass as well as massive galaxies, hence the 
problem of warm ISM content in galaxies will also be addressed. 

The paper is structured as follows: in Sect. 2 we introduce the 
sample and in Sect. 3 we discuss the methods of the data analysis. 
In Sect. 4 we overview various contaminating factors polluting un- 
resolved emission. The properties of unresolved emission in late- 
and in early-type galaxies are discussed in Sect. 5. The anomalous 



emission in the hard band identified in four massive galaxies in 
Virgo cluster is discussed in Sect. 6. We conclude in Sect. 7. 



2 SAMPLE SELECTION 

Since we focus on the unresolved emission, we require at least 
moderately deep observations allowing us to detect and remove 
majority of bright X-ray binaries. From the Chandra archive we 
selected full-size early-type galaxies having source detection sen- 
sitivities better than 5 x 10'^ erg s"'. However, we excluded bright 
cD ellipticals from our sample, such as M87 or NGC1399, further- 
more we also excluded the peculiar radio galaxy, NGC5128. 

To avoid projection effects we only selected face-on spi- 
ral galaxies with point source detection sensitivity better than 
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10^' erg s"' and we also included three well observed irregular 
galaxies. Note, that in late-type galaxies both HMXBs and LMXBs 
are present, hence the overall contribution of X-ray binaries is 
higher. Therefore a better source detection sensitivity is required to 
avoid substantial pollution of the unresolved emission with X-ray 
binaries - for detailed discussion see Sect. |4.2| 

Our sample consists of 26 galaxies, which includes 1 1 late- 
type and 15 early-type systems. Major properties of the selected 
galaxies are listed in Table [T] 



3 DATA REDUCTION 
3.1 Chandra 

We analyzed all available Chandra observations for the se- 
lected galaxies, which had exposure times longer than 2 ks. In total, 
this yielded a 4.6 Ms data. The data analysis was performed with 
standard CIAO software package tools (CIAO 4.3; CALDB 4.4.2). 

For each observation the flare contaminated time intervals 
were filtered, excluding those intervals where the count rate devi- 
ated more than 20% from the mean value. If a galaxy was observed 
by multiple pointings, we projected and merged each observations 
to the coordinate system of the observation with the longest expo- 
sure time. The combined observed and filtered exposure times for 
each galaxy are given in Table [T] To maximize the source detec- 
tion sensitivity, we ran the CIAO wavdetect tool on the unfiltered 
combined data set in the 0.5 - 8 keV band energy range. As our 
particular goal was to study the unresolved X-ray emission, which 
demands careful removal of point sources, several parameters of 
the WAVDETECT tool Were changed from their default values. The 
scales on which we searched for sources were the V2 series from 
I.O to 8.0. The value of eenergy was set to 0.5, which parameter 
describes the percentage of PSF energy to be encircled. We also 
changed the values of maxiter to 10 and iterstop to 0.00001. The 
value of sigma parameter was set to 4, thereby minimizing the con- 
tribution of residual counts from points sources to the difi'use emis- 
sion. We thus obtained relatively large source cells, which include 
a large fraction of source cells. The resulting list of sources were 
used to mask out point sources for further analysis. 

To estimate and subtract the instrumental and sky background 
components we used a combination of several regions away from 
the galaxies for all galaxies except for M31, M51, M74, M81, M83. 
The angular size of these galaxies is comparable or larger than 
the extent of the combined Chandra image, hence a direct back- 
ground subtraction is not possible. Therefore we used the ACIS 
"blank- sky" files ^http://cxc.harvard.edu/contrib/maxim/acisbg / 1 to 
estimate the background. As the instrumental background compo- 
nents of Chandra vary with time, we renormalized the background 
counts using the 10- 12 keV count rates. The computed source de- 
tection senitivities are listed in Table[T] and they refer to the 0.5 - 8 
keV energy range assuming Galactic column density and a power 
law model with slope of F = 1.56, which is typical for LMXBs 
llrwin et al. 2003). 

To obtain X-ray luminosities of the unresolved emission in 
the sample galaxies we extracted spectra using the CIAO specex- 
TRACT tool. For each galaxy we used the extraction regions given 
in Table [T] The X-ray energy spectra were fit with a two compo- 
nent model, consisting of an optically-thin thermal plasma emis- 
sion model (Mekal model in Xsp ec) with abundances fixed at Solar 
value jAnders & Grevess3ll989 l) and a power law model. The col- 
umn density was fixed at the Galactic value for early-type galaxies. 



whereas it was fixed at Nh = 10^' cm"^ for late-type galaxies to 
approximately account for the additional intrinsic absorbtion. All 
fits were performed in the 0.4 - 7 keV energy range. Due to the 
relatively low number of counts in NGC3377 and NGC3585, the 
temperature of the Mekal model was fixed at kT = 0.40 keV. The 
0.5 - 2 and 2-10 keV band X-ray luminosities of the unresolved 
emission are listed in Table [2] The meaning and significance of 
these values are described and discussed throughout the rest of the 
paper. 

3.2 Infrared data 

We used the 2MASS K-band data jjarrett et ai]|2003l) to trace 
the stellar light in all galaxies but M31. In case of M3I the K-band 
image, provided by the 2MASS archive, suffers from a background 
subtraction problem and is not suitable for the analy sis. Instead, 
we rely on the 3.6 yum Infrared Array Camera (IRAC) l lFazio et al] 
,2004) data of Spitzer Space Telescope. To estimate the background 
level of these images, we used nearby regions off the galaxy. To 
facilitate the comparison with other galaxies, we converted the 
Spitzer 3.6 /jm counts to 2MASS counts. The conversion factor is 
Ci.i ^m-ilC K-band ~ 10.4, that was obtained in the central region of 
M31 where the role of background is negligible. 

To measure star-formation rates (SFRs) in late-type galaxies 
we u sed the 70 ^m data of Multiband Imaging Photometer (MIPS) 
(Rieke et alj|2004l) onboard of Spitzer Space Telescope. The far- 
infrared data was con verted into total infrared luminosity based on 
iBavouzet etal] ( |2008|) . whi ch va l ue wa s used to compute the corre- 
sponding SFR according to lBelll ^20031) . 



4 SOURCES OF CONTAMINATION OF UNRESOLVED 
EMISSION 

In the soft and hard X-ray bands a number of contaminating 
factors pollute the unresolved emission, whose contribution needs 
to be subtracted. In this section we overview these factors, and their 
effect on the derived X-ray luminosities. 

4.1 Residual counts from resolved compact sources 

Although we use relatively large source cells to exclude the 
resolved sources, a certain fraction of source counts falls outside 
these regions. Their contribution must be removed when studying 
unresolved X-ray emission. 

To calculate the residual emission from resolved sources we 
employed the following procedure. For each source we extracted 
the point spread function using CIAO mkpsf tool and computed 
the fraction of counts falling outside the source cell. For most of 
the sources this fraction is ~2%, if it was larger, the source cell 
was enlarged accordingly. Few galaxies in our sample, for ex- 
ample NGC4278 or M81, hosts extremely bright central sources 
(Lx > lO"*" erg s"'), which could significantly influence the ob- 
served unresolved X-ray emission. In these cases, the source cells 
were enlarged to contain >99% of the source counts. The resid- 
ual counts from all sources were summed and subtracted from the 
unresolved emission. Their contribution varied between 1-25%. 

4.2 Unresolved X-ray binaries 

For the purpose of this study we consider all X-ray binaries, 
including um'esolved ones, as a source of contamination. Resolved 
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Table 2. Observed X-ray luminosities and LxlLf; ratios of the sample galaxies in the 0.5 - 2 and in the 2 - 10 keV energy range. 



Name 


kT 


^0.5- 


2ki;V 


-f'2-IOkeV 


io.5-2kcV,XB,.sub 


i2-10keV,XB,sub 


i'2-10ki;V,XB,ISM,sub 


-^■0.5-2keV,XB,sub / I^K 


-^■2-10kcV,XB,ISM,sub/i'A: 




(keV) 


(erg 


s-') 


(erg s"') 


(erg 




(erg s"') 


(erg 




(erg s ' 


LkIo) 


(ergs ' 


LkIo) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


M51 


0.32 


5.3 X 


103' 


1.9 X 103' 


5.0 X 


1039 


1.5 X 1039 


1.4 X 


1039 


7.9 X 


1028 


2.3 X 


1028 


M74 


0.24 


4.7 X 


1038 


3.5 X 1038 


4.1 X 


1038 


2.3 X 1038 


2.3 X 


1038 


2.5 X 


1028 


1.5 X 


1028 


M81 


0.32 


3.5 X 


1038 


2.6 X 1038 


3.4 X 


1038 


2.5 X 1038 


2.5 X 


1038 


7.1 X 


1027 


5.2 X 


1027 


M83 


0.34 


3.4 X 


103' 


7.6 X 1038 


3.3 X 


1039 


5.9 X 1038 


5.8 X 


1038 


8.4 X 


1028 


1.5 X 


1028 


M94 


0.47 


1.2 X 


103' 


3.9 X 1038 


1.2 X 


1039 


3.2 X 1038 


3.2 X 


1038 


3.6 X 


1028 


9.6 X 


1027 


M95 


0.24 


7.4 X 


1038 


4.0 X 1038 


6.7 X 


1038 


2.8 X 1038 


2.8 X 


1038 


2.0 X 


1028 


8.5 X 


1027 


MlOl 


0.22 


7.3 X 


1038 


2.6 X 1038 


7.2 X 


1038 


2.5 X 1038 


2.5 X 


1038 


2.1 X 


1028 


7.4 X 


1027 




0.24 


7.5 X 


1037 


5.1 X 1037 


7.1 X 


1037 


4.4 X 1037 


4.4 X 


1037 


1.4 X 


1028 


8.8 X 


1027 


NGC3077 


0.31 


7.8 X 


1037 


6.5 X 1037 


6.2 X 


1037 


4.9 X 1037 


4.9 X 


1037 


2.4 X 


1028 


1.9 X 


1028 


NGr47 1 4 


0.18 


7.5 X 


1037 


1.7 X 1037 


5.1 X 


1037 


1.1 X 1037 


1.1 X 


1037 


7.9 X 


1028 


1.7 X 


1028 


NGC4449 


0.20 


5.5 X 


1037 


1.6 X 1038 


5.3 X 


1038 


1.4 X 1038 


1.4 X 


1038 


1.6 X 


1029 


4.1 X 


1028 


M31 bulge 


0.32 


2.9 X 


1038 


1.8 X 1038 


2.9 X 


1038 


1.6 X 1038 


1.6 X 


1038 


7.8 X 


1027 


4.3 X 


1027 


M32 


0.60 


3.0 X 


1036 


3.3 X 103* 


3.0 X 


1036 


3.3 X 103* 


3.3 X 


103'5 


3.5 X 


1027 


3.9 X 


1027 


M49 


0.82 


1.1 X 


lo-" 


2.3 X Iff*" 


1.1 X 


lo-" 


1.8 X 10'*" 


1.1 X 


1040 


3.6 X 


1029 


3.8 X 


1028 


M60 


0.73 


6.9 X 


1040 


1.1 X 10**" 


6.8 X 


1040 


8.8 X 1039 


5.9 X 


1039 


3.1 X 


1029 


2.7 X 


1028 


M84 


0.61 


3.4 X 


1 n40 
10 


5.5 X 1039 


3.4 X 


1 n40 
10 


3.3 X 1039 


2.4 X 


10 


2.1 X 


10 


1.5 X 


10 


M89 


0.54 


1.7 X 


1040 


2.6 X 1039 


1.6 X 


1040 


5.3 X 1038 


2.3 X 


1038 


2.3 X 


1029 


3.2 X 


1027 


M105 


0.64 


1.6 X 


1038 


1.1 X 1038 


1.5 X 


1038 


8.6 X 1037 


8.6 X 


1037 


3.7 X 


1027 


2.1 X 


1027 


NGC1291 


0.31 


1.3 X 


1039 


4.5 X 1038 


1.3 X 


1039 


2.6 X 1038 


2.6 X 


1038 


2.0 X 


1028 


4.1 X 


1027 


NGC3377 


0.4 


1.8 X 


1038 


3.3 X 1038 


6.0 X 


1037 


3.8 X 1037 


3.8 X 


1037 


3.0 X 


1027 


1.9 X 


1027 


NGC3585 


0.4 


3.1 X 


1039 


2.4 X 1039 


2.3 X 


1039 


5.4 X 1038 


5.4 X 


1038 


1.5 X 


1028 


3.6 X 


1027 


NGC4278 


0.47 


9.1 X 


1038 


3.7 X 1038 


8.3 X 


1039 


1.8 X 1038 


1.7 X 


1038 


1.5 X 


1028 


3.1 X 


1027 


NGC4365 


0.49 


3.5 X 


1039 


2.0 X 1039 


2.8 X 


1039 


3.4 X 1038 


3.1 X 


1038 


2.6 X 


1028 


2.8 X 


1027 


NGC4526 


0.33 


3.0 X 


1039 


1.9 X 1039 


2.4 X 


1039 


3.4 X 1038 


3.3 X 


1038 


2.5 X 


1028 


3.4 X 


1027 


NGC4636 


0.59 


1.2 X 


lo-" 


6.9 X 1039 


1.2 X 


lo-" 


5.1 X 1039 


2.5 X 


1039 


1.5 X 


1030 


3.1 X 


1028 


NGC4697 


0.33 


1.0 X 


1039 


2.9 X 1038 


9.5 X 


1038 


8.6 X 1037 


8.3 X 


1037 


1.9 X 


1028 


1.6 X 


1027 



Note. Columns are as follows. (1) Best-fit temperature of the soft component described by Mekal model in Xspec. In case of NGC3377 and NGC3585 the 
temperature was fixed at the given value. (2) and (3) Total observed unresolved X-ray luminosity in the 0.5 - 2 keV and 2-10 keV band, respectively. The 

CXB contribution and the contribution of source counts falHng outside the source cells is subtracted. (4) and (5) X-ray luminosity after the emission from the 
population of unresolved X-ray binaries are subtracted in the 0.5 - 2 keV and 2-10 keV band, respectively. (6) Unresolved X-ray luminosity in the 2-10 

keV band after the emission from warm ISM is subtracted. (7) and (8) Final, contamination subtracted, LxILk ratios in the 0.5 - 2 keV and 2-10 keV band, 

computed from the X-ray luminosities listed in columns (4) and (6) respectively. 



point sources account for the bulk of the emission from X-ray bina- 
ries and their removal is rather straightforward. They were masked 
out and excluded from the further analysis as a part of the data 
preparation procedure. The issue of the "spill-over" counts is ad- 
dressed in the previous section. Unresolved X-ray binaries, on the 
contrary, cannot be removed on the source-by-source basis. As 
they make a notable contribution to unresolved emission in some 
galaxies, their contribution has to be removed statistically, based 
on the knowledge of their luminosity distributions and Lx/Lk and 
Lx/SFR ratios. Therefore, we distinguish high-mass and low-mass 
X-ray binaries. The task is further complicated by the fact that 
LxIL^ and L^/SFR ratios for LMXBs and HMXBs are not exactly 
constant and may vary from galaxy to galaxy. 

In early-type galaxies, which host only LMXBs, we circum- 
vent this difficulty by determining the normalization of the X-ray 
luminosity function (XLF) from the number of resolved X-ray bi- 



naries. In this computation we took into account that certain frac- 
tion of detected sources are not X-ray binaries but resolved cos- 
mic X-ray backgroun d (CXB) so urces. We computed their num- 
ber based on M oretti et al. I ( l2003h . and found that in all galaxies 
LMXBs comprise the majority, k80 - 97%, of compact sources. 
We also assumed that the shape of the XLF of LMXBs is described 
by their average luminosity function derived in .Gilfanov. ( ,2004]) . 
To transform it to the 0.5 - 2 keV and 2-10 keV energy bands we 
used the ave rage LMXB spec trum, namely a power law with slope 
of r = 1.56 l llrwin et alj2003h . To compute the luminosity of unre- 
solved binaries, the luminosity functions were integrated down to 
103^ erg s"'. The thus obtained luminosity of unresolved LMXBs 
was subtracted from the observed luminosity of each galaxy (Table 

In most galaxies the contribution of LMXBs in the 0.5 - 2 
keV band does not exceed 20%, except for NGC3377 (67%) and 
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Figure 1. Top panels: X-ray versus K-band luminosity for a set of elliptical, spiral, and irregular ga laxies in the 0.5 - 2 k eV band (left panel) and 2-10 keV 
band (right). The dashed line shows total emissivity of faint sources in the Solar n eighborhood f rom Sazon ov et al The solid line (red) is the same but 

excluding contribution of young stars. In the 0.5 - 2 keV band we recalculated the lSazonov et a l. 1 2006) LxILk ratios from the 0.1 - 2.4 keV value (Sect. l5.lT 
Bottom panels: Same data presented in the form of LxILk ratios. The solid and dashed lines represent the same as in the upper panels. 



NGC3585 (26%). In the 2- 10 keV band unresolved LMXBs play a 
more significant role, their contribution typically varying between 
19 - 50%. But, in six galaxies with low LxIL^ ratios ~ namely 
M89, NGC3377, NGC3585, NGC4365, NGC4526, NGC4697 - a 
large fraction, 70 - 89%, of the total observed emission is due to 
unresolved LMXBs (Table [2}. Therefore the accuracy of LxlLf^ 
ratios of these galaxies may be somewhat compromised by residual 
emission from unresolved LMXBs, because of the limited accuracy 
of the employed procedure. 

In late-type galaxies both LMXBs and HMXBs are present, 
which - by means of X-ray observations ~ cannot be distinguished. 
For this reason, their X-ray luminosity functions cannot be renor- 
malized individually. We therefore emplyed the f ollowing proce 



J ^—\. 

dure. We used the luminosity func t ions of LMXBs lGilfanovl2004 ) 
and HMXBs jGrimm et a"lll2003l : [shTykovskiv & Gilfanovll200! 



iMineo et al.ll2011 ) with their average normalizations to predict the 
numbers of low-mass and high-mass X-ray binaries above the de- 
tection threshold based on the K-band luminosity and SFR of each 
galaxy, respectively. We then used the ratio of the detected number 
of compact sources (minus predicted CXB contribution) to the pre- 
dicted number to compute the correction factor. This factor was 
used to correct normalizations of LMXB and HMXB XLFs for 
each galaxy. The number of CXB sou rces was predicted based on 
source counts of iMoretti et al. I i2003h . The thus corrected XLFs 
were used to compute the luminosities of um'esolved X-ray bina- 



ries. In converting the luminosities of LMXBs to the soft and hard 
bands we used spectral parameters specified above. For HMXBs 
we assumed a power law model wi th slope of F = 2, corresponding 
to the median value determined bv lSwartz et al.l 12004 ). As before, 
for HMXBs we used Nm = 10^' cm~^, in order to approximately 
account for the intrinsic absorption present in late-type galaxies. 

Obviously, the method of correcting the XLF normalizations 
for HMXBs and LMXBs by the same factor is not accurate. How- 
ever, no better correction is possible, as Chandra data does not al- 
low an unambiguous separation of the two types of X-ray binaries 
in external galaxies. On the other hand, thanks to the good source 
detection sensitivity demanded in the selection of late-type galax- 
ies, the contribution of unresolved binaries is rather small and this 
inaccuracy does not significantly compromise the obtained LxILk 
ratios. Indeed, in the 0.5-2 keV band the contribution of unresolved 
X-ray binaries does not exceed < 15%. It is somewhat larger in the 
2-10 keV energy range, but still remains below ~ 35% in all 
galaxies (Table|2ll. 



4.3 Contribution of warm ISM to the hard band 

Although the warm ISM has a temperature in the kT = 0.2 - 
0.8 keV range and bulk of its emission comes below 2 keV, in some 
gas rich galaxies it may contribute to the hard band. To be precise. 
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this component should not be considered as a contaminating factor. 
However, as one of our aims is to study the 2-10 keV emission 
from faint unresolved sources, it is desired to subtract the contribu- 
tion of warm ISM from the observed hard band emission. 

This calculation is rather straightforward. For the Mekal 
model with the temperature fixed at the best-fit value (Table |2j we 
compute the (0.5 - 2 keV)/(2 - 10 keV) hardness ratio and use it to 
determine the luminosity of the warm ISM in the 2-10 keV band. 
Its contribution is virtually negligible for the majority of galaxies 
except for M49, M60, M84, M89, and NGC4636 where it can ac- 
count for up to ~ 10 - 30% of the observed hard band emission 
(TableO. Note, that the employed procedure slightly overestimates 
the contribution of warm ISM to the hard band, because we used 
the total 0.5 - 2 keV luminosity in this calculation. However, in 
these five galaxies the soft band emission is largely dominated by 
the ISM contribution, therefore this inaccuracy does not influence 
the final result in any significant way. 



5 X-RAY TO K-BAND LUMINOSITY RATIOS 

5.1 Observed LxILk ratios 

The final, contamination subtracted, Lx I Lie ratios are listed in 
the last two columns of Table |2] for the soft and hard bands, re- 
spectively. Additionally, the contamination subtracted X-ray lumi- 
nosities as a function of K-band luminosities are d epicted in Fig. 
T] wh ere the Solar neighborhood Lx - Lx relation dSazonov et al.l 
2006 ') is also shown. The dashed line shows the correlation for all 
types of observed sources, including CVs, ABs, and young stars, 
whereas in drawing th e solid line we excluded the contribution of 
young stars. Note that lSazonov et aL I j2006l) computed Lx/Lk ra- 
tios in the 0.1 - 2.4 keV band, hence we converted their values into 
the 0.5 - 2 keV energy range assuming a power law model with 
slope of r = 2. 

Ideally, gas-free early-type galaxies are expected to lie along 
the solid line in Fig. [T] Late-type systems, on the contrary, are not 
neccessarily expecte d to lie along the dash ed line for two reasons. 
On the one hand, the lSazonov et al.l ( |2006[) relation was calibrated 
for the Solar neighborhood; on the other, the population of young 
stars is expected to correlate with the star-formation rate rather than 
the stellar mass (Sect.[53). 

Although there is an obvious correlation between the X-ray 
luminosity and stellar mass of the galaxy, there is large scatter in 
the Lx /Lk ratio in both energy bands and in galaxies of all morpho- 
logical types. The observed scatter is smaller for low-mass galax- 
ies, but more massive systems (Lk ^ 5 x 10'" Lk.o) tend to have 
10 - 100 times higher X-ray luminosities than expected based on 
their K-band luminosities. 



5.2 Early-type galaxies 

It is well known that there is a broad correlation between the 
stellar mass of early-type galaxies and their soft X-ray luminosity, 
which relation bears a lar ge dispersion and is non -linear, especially 
at the high-mass end (e.g. io'Sullivan et al.ll200ll) . Our results fit in 
this picture. The four least massive galaxies with Lx < 4x 10'" Lk,g 
have Lx/Lk ratios consistent or even smaller than the Milky Way 
value, indicating their almost zero gas content. At higher galaxy 
mass the Lx / Lk increases, exceeding the Milky Way value by factor 
of 3 - 250 and reaching values characteristic of massive gas-rich 
ellipticals. 



10^ 



10' 



10' 



Early-type )K 

























Age, Gyr 



10 



12 



Figure 2. Soft band Lx/Lk ratios as a function of the age for gas-rich early- 
type galaxies with K-band luminosity exceeding 5 X lO'" Lq. Note, that 
NGC1291 is not included since no reference for its age has been found. 



The main reason for the relatively gas-free nature of low-mass 
ellipticals is their shallow potential well, that permits the formation 
of SN la driven galactic-scale outfl ows, which are capable of re- 
moving major fraction of the ISM ( David et al ] |2006l : iLi & Wan^ 
2007; Bogdan & Gilfanov 2008). Massive systems have a deeper 
potential well, and the energy supplied by SNe la is not sufficient 
to launch an outflow, hence the ISM will be accumulated in the 
galaxy. 

Calculations suggest that an important factor influencing the 
gas content of a galaxy is its age (e.g. O'Sullivan & Ponman 2004i), 
younger galaxies having less gas. Our results appear to support this 
conclusion. Indeed, three rather massive galaxies at Lk ~ 10" Lq, 
NGC3585, NGC4365, NGC4526, host fairly small amount of ISM . 
Their stellar ag es are young: 3. 1 Gyrs jTerlevich & Forbes" 2002h. 
3.6 G yrs (De nicolo et al.ll2005b . and -1.7 Gyrs (Gallagher et al] 
l2008h . respectively. On the other hand, all gas rich massive galax- 
ies in our sample have ages of ~I0 Gyrs (e.g. iTerlevich & Forbes! 
l2002h . A possible existence of some correlation between Lx/Lk ra- 
tio of massive galaxies with their age is further illustrated by Fig. 
|2l where we plot all early-type galaxies with K-band luminosity 
exceeding 5 x 10'" Lq, except for NGC1291 for which no age ref- 
erence has been found. The obtained correlation is not very tight, 
which is not surprising as a number of other parameters may play 
a role, such as the environment (see Sect. ^ and the merger his- 
tory of the galaxy. As another caveat, we mention that the stellar 
age of galaxies used for this plot were determined for central part 
of galaxies, typically inside r(,/8 (r^ is the effective radius, contain- 
ing 50% of the stellar light), which corresponds to a few arcsec for 
the galaxies in our sample. Therefore, the listed age measurements 
may be inaccurate for galaxies with strong age gradients. 

Three low-mass galaxies, M32, M105, NGC3377, have Lx/Lk 
ratios in the soft band at the « (3-4)x 10"^ erg s"' Lj^'^ level, which 
is by a fac tor of ~2 lower than t he Solar neighborhood value deter- 
mined by ISazonov et alj ( |2006|) . Two of these galaxies, M32 and 
NGC3377 are presumabl y gas-free, whereas M l 05 may contain a 
small amount of hot gas dTrinchieri et al.ll2008h . Note, that in case 
of NGC3377 the relatively poor source detection sensitivity may 
somewhat influence the obtained Lx/Lk ratio (Sect. |42] |. The ob- 
served factor of ~2 difference may be a result of some residual con- 
tamination by young population in the Solar neighborhood value, or 
caused by inaccurate spectral band conversion. In principle, it may 
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Figure 3. X-ray luminosity in the 2-10 keV band as a function of K- 
band luminosity for early-type galaxies. The contribution of warm ISM and 
unresolved LMXBs is subtracted. The galaxies located in Virgo cluster are 
marked with the large (red) symbols. 




Figure 4. Unresolved X-ray luminosity in the 2- 10 keV band as a function 
of the SFR for spiral and irregular galaxies. The solid line shows the relation 
Lx/SFR= 1.7 X 10'^ (erg/s)/(Mo/yr) which is about 7% of the total value 
associated with HMXBs. 



also be a consequence of the galaxy-to-galaxy variations, but this 
possibility seems less likely, in the view of the hard band LxlLf^ 
ratios (Fig.lTJ. 

In the 2-10 keV energy range, the majority of early-type 
galaxies (11 out of 15) show fairly uniform LxILk ratios, fluctu- 
ating around the Solar neighborhood value. The remaining four 
galaxies, massive ellipticals in the Virgo galaxy cluster M49, M60, 
M84, and NGC4636, are significantly more luminous in the hard 
band. The statistical uncertainty of the Lx/^a: ratio determination 
for these galaxies is (0.2 - 0.4) x 10^^ erg s"' L^'^, therefore the 
deviation is highly significant from the statistical point of view. 
Excluding these four galaxies, we calculate for the 1 1 early-type 
galaxies an average ratio of Lx/^a: = (3.1 ±0.9) x 10^' erg s"' Lj^'^, 
where the cited error is the rms of the individual values. This num- 
ber is in excellent agreement with the (3.1±0.8)xlO^'' erg s" ' Lj^'^ 
obtained by Sazonov et al. ( 2006) for the Solar neighborhood. In 
iBogdan &^Gilfanovl ^2010 ) we demonstrated that the radial surface 
brightness profiles of gas-poor ellipticals in the soft X-ray band fol- 
low the stellar light distribution, and that spectral charachteristics 
of these galaxies agree well with each other. These facts suggest 
that the bulk of the hard band X-ray emission is due to faint un- 
resolved stellar X-ray sources associated with old populations - 
predominantly active binaries and cataclysmic variables. 

Interestingly, the four galaxies having enhanced LxILk ratios 
in the hard band, belong to Virgo cluster of galaxies (Fig. [3}. The 
possible origin of the excess hard X-ray emission in these galaxies 
is discussed in Sect. [6] 

5.3 Late-type galaxies 

Late-type galaxies tend to have systematically larger X-ray 
luminosity in both bands, than ellipticals of the same mass (Fig. 
[Hi. For low-mass and medium-mass galaxies this is not too sur- 
prising, as late-type galaxies host additional emission components 
- gas and young compact sources, associated with star formation. 
A similar comparison with massive ellipticals cannot be made be- 
cause the mass of late-type galaxies in our sample is limited by 
Lx ~ T X 10'" Lq. The Lx/Lx ratios of late-type galaxies exhibit 
a large scatter, not only in the soft band but also in the hard band. 
This can be caused by three reasons: (i) similarly to ellipticals, the 



gas content of late-type galaxies may vary from galaxy to galaxy; 
(ii) the population of young stars and YSOs must be determined 
by the star-formation rate rather than by the mass of the galaxy, re- 
sulting in the SFR-dependent Lx/Lx ratios; (iii) the Lx/Lx ratio of 
the populations of CVs may depend on the average age of the stel- 
lar population, which may be rather young and varying in irregular 
galaxies and disks of spiral galaxies. 

The X-ray energy spectra of all late-type galaxies in our 
sample shows an enhance d soft component as demonstrated in 
iBogdan & Gilfanovl ( I20T J). This component can be described by 
an optically-thin thermal plasma emission spectrum with the tem- 
perature in the range of kT ~ 0.2-0.5 keV (Table[2j, suggesting the 
gaseous nature of the emission. It is plausible to conclude that all 
late- type galaxies contain at least moderate amount of ionized gas. 
Variations in the amount of hot gas may account partly or entirely 
for the observed scatter in the Lx/Lx ratios in the 0.5-2 keV en- 
ergy range. The contribution of this component is negligible above 
~2 keV. 

In the 2-10 keV band the high Lx/Lx ratios are caused 
by the contribution of young stars a nd YSOs which are k nown to 
be sources of hard X-ray emission l lKovama et alj[l99^ . As the 
population of these sources is determined by the star-formation 
rate, a correlation of the hard band luminosity with the SFR is 
expected. In Fig. |4] we depict the 2-10 keV band luminos- 
ity of the young population as a function of the SFR. In order 
to remove the contribution of the X-ray sources associated with 
the old population, we subtracted the luminosity corresponding to 
Lx/Lx = 3.1 X 10^' ergs"' L|^'g from the observed values. Fig. 
[4] shows that a correlation with SFR exists indeed, albeit with 
some scatter. For our sample we calculated the average ratio of 
Lx/SFR = (1.7±0.9)x 10'** (erg/s)/(Mo/yr), where the error refers 
to the rms of the observed values. This relation is shown in Fig.|4]by 
the solid line. We also performed a two parameter fit to the data in 
the form Lx = axLx + bx SFR and determined the following best- 
fit values for the scale factors: a = (2.1 ± 2.8) x 10^' erg s"' L^'^^ 
and b = (2.5 ± 0.6) x 10^** (erg/s)/(MG/yr). While the scale factor 
b is consistent with the Lx/SFR value derived above, the uncer- 
taintiy of the factor a is too large to make any meaningful con- 
clusions regarding the average Lx/Lx ratio in late-type galaxies. 
Note, that the obtained L^/SFR ratio for unresolved young stars 
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and YSOs correspon ds to ~7% of the total X-ray luminosity asso- 
ciated with HMXBs jGrimm et alJliooil : Ishtvkovskiv & Gilfanovl 
l2005l:lMineo et^laOlU) ! 



6 EXCESS EMISSION IN THE HARD BAND IN 
MASSIVE VIRGO ELLIPTICALS 

A major unexpected result of our study is the detection of sig- 
nificantly enhanced LxIL^ ratios in the hard band in four massive 
ellipticals (M49, M60, M84, NGC4636). 

Although the soft band LxILk ratio is known to vary signif- 
icantly, due to varying hot gas content of elliptical galaxies, the 
typical gas temperatures even for the most massive galaxies in our 
sample do not exceed a 0.8 keV, therefore the gas contribution to 
the 2-10 keV luminosity is not dominant. The 2-10 keV emis- 
sion from elliptical galaxies is believed to be determined by faint 
unresolved stellar sources, which numbers and total luminosities 
are expected to scale rather uniformly with the mass or K-band 
luminosity of the host galaxy. Indeed, the majority of ellipticals in 
our sample obey this expectation, following the Lx oc Lk law with a 
rather small scatter of ~30%, in contrast to the soft band luminosity. 
This behavior breaks down in the case of the four above mentioned 
galaxies, having approximately 5-12 times larger LxILk ratios in 
the hard band. 



6.1 Common properties 

Main parameters of the Virgo galaxies are summarized in Ta- 
ble[3] listing their LxIL^ ratio, K-band luminosity, age, stellar ve- 
locity dispersion, radial velocity with respect to M87, offset from 
the center of M87, and globular cluster (GC) specific frequency. 
We use M87 to characterize the position and velocity of the galaxy 
with respect to the intracluster medium. From the table, one may 
conclude that galaxies with anomalous Lx/Lk ratios have several 
common properties: (i) are members of the Virgo cluster of galax- 
ies; (ii) are massive, their total K-band luminosities are in the range 
of (1.2 - 3.8) X 10" Lk,o; (iii) have relatively small velocity with 
respect to the intracluster gas Ai\ < 300 - 40 km/s; (iv) are old, 
with the st ellar age in the 8.5- 1 1.8 G yrs range jTerlevich & Forbe3 
l2002l : [sinchez-Blazauez et al.ll2006h . 

Based on our, admittedly limited, sample, it appears that each 
of these properties is essential. Indeed, the following galaxies, 
missing one of the above properties, have perfectly normal Lx/Lk 
ratios: (i) a similarly massive and old Virgo galaxy M89, but having 
a large radial velocity with respect to M87, Ai',- ^ 967 km/s; (ii) a 
massive field galaxy NGC3585 with a total mass of 1.5 x 10" Mq 
and relatively young age; (iii) two massive, L^- ~ 2 x 10" Lq, 
but apparently young Virgo galaxies NGC4365 and NG C4526 with 
the age of 3.6 and 1 .7 Gyrs, respectively ( fPenicolo et al.l [20051 ; 
iGallagher et ai]|2008l) . one of them having small (Av^ ~ 64 km/s) 
and the other large (Avr ~ 859 km/s) velocity. Interestingly, there 
seems to be no dependence on the radial distance from the center of 
Virgo cluster, which may be explained by the somewhat irregular 
spatial distribution of its intracluster gas. 

As one can see from the Table[3] the four galaxies with anoma- 
lously high Lx/Lx ratios have the highest values of the globu- 
lar clusters specific frequency. Therefore their LMXB populations 
have a larger contribution of X-ray binaries formed dynamically in 
globular clusters, which could, hypothetically, result in enhanced 
X-ray luminosity in the hard band, due to unaccounted contribu- 
tion of faint sources formed in globular clusters. This is further dis- 



cussed in the Section [6.4.2l ( and it is shown that GC sources cannot 
explain the anomalous hard emission. 

Because of the limited size of our sample, it is unclear, 
whether the observations presented above reflect true properties or 
are a result of a chance coincidence. With this in mind, we proceed 
with the discussion of various mechanisms which may be respon- 
sible for the excess hard emission. 

6.2 Energy spectra 

In the left panel of Fig. |5] we show the energy spectra of all 
studied galaxies combined in three groups: 1) the galaxies showing 
anomalous Lx/Lk ratios ~ M49, M60, M84, and NGC4636, 2) all 
remaining early-type galaxies and 3) all late-type galaxies. To fa- 
cilitate comparison, the spectra are normalized to L^: = 10" Lk,o 
and are projected to the distance of 16 Mpc. Firstly, Fig. |5] shows 
that at every energy the spectra of the "anomalous" ellipticals is 
higher by 1 - 2 orders of magnitude. Secondly, it is obvious that the 
enhanced hard Lx/Lk ratios are not a result of the soft component 
"leaking" into the 2-10 keV band. Instead, it requires an additional 
hard component in the spectrum. Such a component, if of thermal 
origin, should have a temperature of about 2 keV, as illustrated by 
the right-hand panel in Fig[5] 

We fitted the spectrum with a model consisting of two Mekal 
models and a power law component. The latter accounts for the 
population of unresolved LMXBs and has the slope fixed at F = 
1.56, and the normalization fixed at the value corresponding to the 
estimated luminosity from unresolved LMXBs. The best-fit tem- 
peratures of the thermal components are kTi = 0.74 ± 0.01 keV 
and kT2 = 1.75 + 0.12 keV. The hotter component can be replaced 
by a power law with the photon index of F = 2.09 ± 0.04, giving 
similar fit quality. In both models, the fit is statistically unaccept- 
able, mainly because of the poor approximation at low energies, 
indicating that the spectrum of the soft component is more com- 
plex. 

6.3 Spatial distribution of hard emission 

The surface brightness profiles of unresolved emission are pre- 
sented in Fig|6] In the left four panels we show for each galaxy the 
soft and hard band profiles along with the K-band light distribu- 
tion. The contribution of the warm ISM emission and unresolved 
LMXBs are subtracted from the hard band profiles. The emission 
from the former was computed from the soft band profiles using 
the best fit temperature of the soft component from Table [2] The 
emission from unresolved LMXBs was calculated from the K-band 
light distribution, the normalizations were determined as described 
in Sect. 14.21 The amplitude of these two main contaminating fac- 
tors is investigated in the right hand panels, where we show the 
total hard band emission, before subtracting ISM and LMXB con- 
tributions (but instrumental and cosmic background removed) and 
the distributions of hard band emission from ISM and unresolved 
LMXBs. 

As one can see from the right hand panels, the hard band emis- 
sion exceeds by a factor of about 5-10 the level of contamination 
(except for M84 where unresolved LMXBs contribute ~l/3 - 1/2 
of the total emission in the hard band). This is much larger than the 
possible uncertainty in their determination. This confirms that the 
observed anomalous hard X-ray emission is not the result of these 
contaminating factors. We remind that the contribution of ISM and 
unresolved LMXBs is subtracted from the hard band luminosities 
plotted in Fig. [T] and Fig.[3] 
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Figure 5. Left: Combined energy spectra of galaxies divided into three groups: 1) the galaxies with anomalous hard Lx/Lk ratios - M49, M60, M84, and 
NGC4636; they are referred to as "Massive ellipticals" in the plot legend; 2) all remaining early-type galaxies, named "Normal ellipticals" and 3) all late-type 
galaxies. The spectra are normalized to the K-band luminosity of lO" Lk.o ™d are projected to a distance of 16 Mpc. The spectra are fluxed using the 
telescope efficiency for a power law spectrum with photon index of 2.09. Right: Optically-thin thermal plasma emission (Mekal) model spectra with three 
different temperatures: kT = 0.3 keV (solid red line), kT = 0.7 keV (dotted black), and kT = 2 keV (dashed green), solar abundance. 



Table 3. Major properties of the Virgo galaxies in our sample. 



Name 


^2- 101icV,XB.ISM,.sub / Lk 




Age 


CTc 


Avr 




Sn 




(erg S-' L^|g) 


(Lk,0) 


(Gyrs) 


(km s-') 


(km s-') 


(') 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


M49 


3.8 X 10^* 


3.8x lO" 


8.5" 


293.8 ± 2.8 


310 


264 


5.40 ± 0.57'' 


M60 


2.7 X 10^** 


2.9 X 10" 


11.0" 


335.3 ± 4.4 


190 


195 


5.16 ± 1.20^ 


M84 


1.5 X 10^** 


2.3 X 10" 


11.8" 


283.3 ±2.8 


247 


89 


5.20 ± 1.45^ 


M89 


3.2 X 10^' 


9.8 X 10'" 


9.6" 


252.6 ± 3.3 


967 


72 


2.82 ± 0.57= 


NGC4365 


2.8 X lO^' 


1.9x 10" 


3.6* 


256.1 ±3.3 


64 


319 


3.86 ±0.71= 


NGC4526 


3.4 X 10^7 


1.5 X 10" 


1.7'' 


263.7 ± 18.9 


859 


286 


1.09 ±0.33= 


NGC4636 


3.1 X 10^8 


1.2 X 10" 


10.3"' 


203.1 ±3.5 


369 


609 


8.9 ± 1.2^ 



Note. Columns are as follo ws. (1) Contamination subt racted Lx/Lk ratios in the 2 - 10 keV energy range. (2) Total K-band luminosity from th e 2MASS 
archive. (3) References are: " iTerlevich & Forbeji2002b - * Denicolo et alj )2005h - = lGallagher etS] j2008l) - '^ ISanchez-Blazquez et al] )2006h . (4) Central 
velocity dispersion from HyperLeda catalog (http://leda.univ-lyonl.fr/ 1. ( 5) Radial velocity rela tive to M87. (6) Dis tance from M87. (7) Globular cluster 

specific frequency. Refererences are: = |Peng et"ai] )2008l) - ^ iDirsch et alj )2005h . 



Although there are some differences between the two distri- 
butions, the soft and hard band emissions have roughly the same 
spatial extent, approximately similar to the spatial extent of the K- 
band emission. 



6.4 Possible origin of the anomalous hard emission 

6.4.1 Accretion of hot intracluster gas 

The best fit temperature of the hard component, kT a 1.7 
keV is close to the temperature of the intraclus ter gas in Virgo 
(e.g. lBohringer et ailll994l : llrwin & Sarazin|[l996l) . suggesting that 
the excess emission may be due to intracluster g as accreted into 
the g ravitational potential of a massive galaxy (e.g. Brigh enti et afl 
1200^) . This interpretation is qualitatively consistent with the four 
abovementioned properties of the galaxies with anomalous hard 
emission. Indeed, (i) such mechanism can work only for galaxies 
located in clusters; (ii) galaxies must be sufficiently massive to con- 



fine the ~2 keV gas; (iii) galaxies should move sufficiently slowly 
with respect to the gas; (iv) some time is needed in order to accrete 
sufficiently massive X-ray halo. 

One of the major facts, of such interpretation would have to 
explain, is the short cooling time of the gas. Given the gas density 
and temperature determined from the spectral fits, its cooling time 
inside 1 kpc does not exceed ~0.3 Gyr in any of the four galaxies. 

Another difficulty of this interpretation is the nearly identical 
radial profiles in the soft and hard bands. Indeed, in the hydrostatic 
equilibrium, the hotter gas component is expected to have broader 
spatial distribution than the cooler one, as the scale heights for gas 
with kT a 1.7 keV and a;0.7 keV are by factors of ~2 - 3 different. 
This, however, is not observed - the soft and hard band emission 
appear to fill the same volume. 

On the other hand, the ISM may have a complex multiphase 
thermal structure, with gas of different temperatures present at 
all radii. Existence of such multiphase ISM i n massive ellip tical 
galaxies has been proposed, for example, by iBuotd ( l2002h and 
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Figure 6. Left: Surface brightness profiles in tlie 0.5 - 2 keV and 2-8 keV bands in four elliptical galaxies with anomalous hard LxIL^ ratios (data points 
with error bars). For each galaxy we also show the K-band light distribution (solid histogram) with arbitrary normalization. The contaminating contribution of 
warm ISM and unresolved LMXBs is subtracted from the 2-8 keV band profiles. Right: Points with error bars show the total 2-8 keV band X-ray surface 
brightness profile, before subtracting the contribution of warm ISM and unresolved LMXBs. The radial distributions of the latter two are shown by dashed and 
dotted histograms. The solid histogram shows the distribution of the K-band light. 
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iBuote et alj i2003h . Based on X-ray spectroscopy of NGC1399 
and NGC5044 they suggested that these galaxies host muhiphase 
ISM. The temperature of the cooler component was found to be 
kT = 0.5 - 0.6 keV, whereas the temperature of the hotter compo- 
nent was consistent with the intracluster medium (ICM) tempera- 
ture, i.e. similar to our results. It remains to be seen, whether ac- 
cretion of hot ICM can explain anomalous hard emission from the 
four massive galaxies. A major advance in testing this hypothesis 
can be achieved by means of high resolution X-ray spectroscopy. 

6.4.2 Faint compact objects 

In principle, it is possible that the hard emission is due to an 
enhanced population of faint stellar sources, such as ABs and CVs. 
This would explain the hard spectrum of the excess emission, con- 
sistent with the power law with the photon index of and the 
fact that hard band profiles approximately follow the distribution 
of the K-band light. Although this assumption cannot be ruled out 
entirely, we stress that the numbers of resolved LMXBs in these 
galaxies are in very good agreement with the prediction based on 
the average scaling relations for LMXBs (Gilfanov 2004). The stel- 
lar ages of the galaxies under consideration are also entirely normal 
for elliptical galaxies. For these reasons, it seems implausible that 
in these galaxies the ratio of A'ab.cv /^lmxb is 5 - 12 times higher 
than in other early-type galaxies from our sample. 

The high specific frequency of globular clusters (GCs) ob- 
served in M49, M60, M84, and NGC4636 may suggest, that the 
excess hard band emission could be due to unaccounted contribu- 
tion of undetected faint sources in GCs. These sources may be un- 
detected, but relatively luminous LMXBs, whose luminosities fall 
just below the Chandra source detection threshold, as well as much 
fainter sources with luminosities in the log(Lx) < 31-33 range 
(CVs, quiescent LMXBs, millisecond pulsars etc). This, however, 
does not seem to be a viable possibility for the following reasons. 
The method of the XLF renormalization is based on the numbers 
of all observed bright LMXBs, irrespective of their association 
with GCs. Therefore the possibly enhanced specific frequency of 
LMXBs in these galaxies, due to the contribution of systems dy- 
namically formed in GCs, is automatically taken into account by 
this procedure and is included in the calculation of the luminosity 
of unresolved LMXBs. In fact, the shapes of the luminosity func- 
tions of field and GC LMXBs are markedly difi^erent: in the latter 
the fraction of fai nt sources (log(Lx) < 37) is ~4 times smaller 
jZhang et alj201 ih . Therefore this procedure may slightly oversub- 
tract the contribution of undetected GC LMXBs. It is also worth to 
note that the specific frequency of resolved X-ray binaries in these 
galaxies is not much different, within a factor of ~1.5, from other 
galaxies in our sample, which fact speaks against significantly en- 
hanced population of (bright) X-ray binaries. 

The possibly enhanced population of faint sources (log(Lx) < 
33) caimot be measured directly in external galaxies but their con- 
tribution can be constrained. From Chandra observations of GCs 
in the Milky Way, the total luminosity of faint sources even in 
those of them with the la rgest encounter rates, e.g. NGC6440 and 
47 Tuc jPoolev et alj|2003 ) does not exceed few x 10^^ erg/s. The 
latter number was obtained by integrating the lumi nosity functions 
of faint sources in NGC6440 jPoolev et al1l2002l) and in 47 Tuc 
( iHeinke et alj2 005') and refers to the 0.5 - 2.5 keV band. The num- 
bers of globular clusters in massive ellipticals is fewx 10^, therefore 
their total luminosity in the soft band does not exceed ~I0^' erg/s. 
The soft-to-hard band ratios for faint globular cluster sources are in 
the ~I - 10 range, but brighter sources tend to have softer spectra 



( iHeinke et alj2005h . Thus an upper limit on the hard band luminos- 
ity is in the ~10'* - 10^' erg/s range. This is about 2-3 orders of 
magnitude smaller than the hard band excess luminosity detected 
in galaxies under consideration, ~I0'''' - lO''" erg/s. 

An even much less likely possibility is that the enhanced 
X-ray emission is due to population of young starts and YSO. 
Using the scaling relation Lx/SFRa; 1.7 x lO'** (erg/s)/(Mo/yr) 
obtained in Sect. 15.31 we conclude that a star-formation rate of 
~ 15-60 Mo/yr is required in order to maintain 2-10 keV lu- 
minosity of Lx = (2 - 1 1) X 10^' erg s"' (Table|2]l. This values are 
appropriate for intensely star-forming g alaxies and are unrealistic 
for elliptic als of « 8.5 - 11.8 Gyrs age jTerlevich & Forbe3l2002l : 
.Sanchez-Blazguez et alj200^ . 



7 CONCLUSION 

We investigated the properties of unresolved X-ray emission 
in a broad sample of nearby early-type and late-type galaxies based 
on archival Chandra data. After removing the contribution of re- 
solved and unresolved X-ray binaries we measured LxILk ratios in 
the 0.5 - 2 keV and 2-10 keV bands and compared them with 
the Solar neighborhood values. Complementing this data with the 
spectral and spatial information we concluded that the unresolved 
X-ray emission originates from, at least, four distinct components. 

(i) The population of faint unresolved sources associated with 
old stellar population. Based on the data for elliptical galaxies we 
obtained a scaling relation for its 2-10 keV luminosity: LxIL^ « 
3.1 X 10"' erg s"' L^'^ with the rms* 0.9 x 10^^ erg s"' L^^, which 
number is in excellent agreement with the value obtained by 
ISazonov et alj ^20061) for the Solar neighborhood. The L^/L^ for 
the soft band and for the late-type galaxies cannot be determined 
unambiguously because of the contribution of ISM and young stel- 
lar sources. For three low-mass gas poor ellipticals we obtained in 
the soft band a;(3 - 4) x 10^^ erg s"' L^'^, which is by a factor of 
~2 smaller than the Solar neighborhood value. 

(ii) In all galaxies, warm ISM with the temperature in the kT ~ 
0.2 - 0.8 keV range is present. In our sample, the gas temperatures 
in late-type galaxies tend to be lower than in early-types galaxies 
of t he same mass. Similar to the results found in previous studies 
(e.g. lO' Sullivan et al fcooil) . the amount and luminosity of the gas 
generally scales with the stellar mass of the host galaxy, albeit with 
the large scatter. The scale factor and the scatter appear to increase 
with the mass of the galaxy. 

(iii) The population of unresolved young stars and YSOs in 
late- type galaxies. The X-ray emission of this component approx- 
imately scales with the SFR of the host galaxy with the average 
Lx/SFR= (1.7 ± 0.9) X (erg/s)/(Mo/yr). 

(iv) The most unexpected result of our study is the detection 
of anomalous emission in the 2-10 keV band from four old and 
massive Virgo ellipticals (M49, M60, M84, NGC4636). We could 
not offer an unambiguous explanation of this emission. A plausible 
interpretation may be that it is a result of accretion of the intraclus- 
ter gas in the gravitational potential of the massive galaxy. In this 
scenario, the two component spectra on one hand and the similar- 
ity of surface brightness profiles on the other, would point at the 
multiphase nature of ISM in these galaxies. High resolution X-ray 
spectroscopy may shed further light on the origin of the anomalous 
emission. 
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